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Abstract. We study the Einstein Yang-Mills Higgs equations in the 50(3) representation on 
a isotropic and homogeneous flat Universe, in the presence of radiation and matter fluids. We 
map the equations of motion into an autonomous dynamical system of first-order differential 
equations and we find the equilibrium points. We show that there is only one stable fixed 
point that corresponds to an accelerated expanding Universe in the future. In the past, 
instead, there is an unstable fixed point that implies a stiff-matter domination. In between, 
we find three other unstable fixed points, corresponding, in chronological order, to radiation 
domination, to matter domination, and, finally, to a transition from decelerated expansion 
to accelerated expansion. We solve the system numerically and we confirm that there are 
smooth trajectories that correctly describe the evolution of the Universe, from a remote past 
dominated by radiation to a remote future dominated by dark energy, passing through a 
matter-dominated phase. 
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1 Introduction 

The quest for a theoretically sound explanation for the current acceleration of the expansion 
of the Universe is one of the major challenges of modern theoretical physics. Most proposed 
models can be roughly divided in two groups according to whether they modify gravity or 
introduce some unusual matter content in the form of one or more scalar fields with appro¬ 
priate potentials (for references see e.g. [1, 2]). In addition to these, there are interesting 
models based on non-local quantum field theory effects [3]. 

The advances in observations have progressively constrained the parameter space of dark 
energy models and, so far, ACDM is still the most favoured one [4]. In this model, dark energy 
is sourced by some vacuum energy that appears as a constant term A in the gravitational 
action and whose effect is, in the present epoch, contrasted by the cold dark matter component 
(CDM). The measured value of A is several orders of magnitude smaller than any realistic 
estimation from quantum field theory (however, see [5] for possible explanations) therefore it 
appears as a very unnatural and fine-tuned explanation for the acceleration of the Universe. 
Future observations (such as the Euclid mission [6]) should be hopefully able to reveal if the 
cosmological constant is truly a constant or it is a slowly varying field of some sort. 

In this paper we assume that dark energy is not a constant vacuum energy but it arises 
as a dynamical effect. Inspired by earlier work by Caldwell et al. [7], and by the multifield 
dynamics studied in the context of Higgs inflation [8], we have already explored in the past 
models of dark energy where gravity is implemented by the dynamics of gauge and Higgs 
fields [9, 10]. In particular, in [10] we have shown that an accelerated phase is possible when 
gravity is minimally coupled to a SU{2) gauge fields and to a Higgs-like complex doublet. 

The Einstein Yang-Mills Higgs (EYMH) equations in a homogeneous and isotropic Uni¬ 
verse were studied in few papers appeared before the discovery of the cosmic acceleration 
[11, 12], so there was no motivation to find a source of dark energy. On the opposite, a lot 
of work has been done in inflationary models driven by gauge fields, eventually coupled to 
the Higgs field [13-15]. Finally, gravity non-minimally coupled to the Higgs field has been 
thoroughly investigated in the context of inflation [16] and compact astrophysical objects [17]. 

In our previous work [10], we did not take in account the coupling between the Higgs 
doublet and the gauge fields. In the present paper, we wish to include this coupling in the 
dynamical system to explore in full generality the physics of the EYMH action on a cosmo¬ 
logical isotropic and homogeneous background. We also choose the representation 50(3), 
instead of SU{2) as we did in [10], mainly to assess whether the late acceleration is some sort 
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of gauge artefact or not. In fact, it is not, the acceleration is real and, as we will shortly see, 
it corresponds to the only stable fixed point of the system of equations. 

One of the conceptual advantages of this model is that it requires degrees of freedom 
whose dynamics is similar to the one of the standard model (SM). Therefore, there is no need 
to modify gravity or to introduce scalar fields with designer potentials that are quite unnatural 
in the realm of particle physics. Here, the standard “Mexican hat” potential, analogous to 
the one of the SM Higgs field, is sufficient to lead to a final accelerated phase. However, as 
we will show, since there are no symmetry breaking effects, the shape of the potential can 
be relaxed to be some generic quartic potential. As already noted in [10], the Mexican hat 
potential does not offer a “slow-roll” phase to the dynamics of the Higgs, it is simply way too 
steep to do that. However, there is another dynamical regime, dubbed “ultra-slow roll” in 
[10], which is related to the kinetic energy of the Higgs phases and can explain the late-time 
domination of dark energy. In fact, the multifield dynamics of the Higgs components modifies 
the Klein-Gordon equation in such a way that it prevents the Higgs field to reach its vacuum 
expectation value in a finite time. Thus, only in the infinite future the potential effectively 
vanish and the Higgs becomes constant. In the meanwhile, the potential behaves as a very 
slowly varying effective cosmological constant equal to the square of the displacement of the 
Higgs field from its vacuum. In the present paper, we will see that the same mechanism is 
at work in the £'0(3) representation, provided we keep account of the coupling between the 
Higgs triplet and the gauge fields. 

The plan of the paper is the following. In the next section we display the action and 
the equations of motion that will be studied as a dynamical system in Sec. HI. In Sec. IV 
we solve numerically the system for realistic initial conditions and we show that dark energy 
can be sourced by the dynamics of the EYMH equations. We conclude in Sec. V with some 
considerations. 


2 EYMH equations in FLRW spacetime 


Let us begin by considering the equations of motion obtained from the Lagrangian 


L = Vi deig\ 

- i(F^$“)(F)^$“) - V($“$“) 

“t“ Ljyi 

, (2.1) 

where 

F;, = d^Al-d,Al + ge-^-AlAl, 


(2.2) 

and 



(2.3) 

is the covariant derivative with the coupling constant g. The term denotes the standard 

Lagrangian of radiation and matter fields in the form of perfect fluids. We choose the repre¬ 
sentation £0(3) so latin indices run over (1,2,3) and summation is understood. As for the 
potential, we choose the standard “Mexican hat” profile 


F = 1 (<1,2 - <1,2)2 ^ 


(2.4) 


where, from now on, <I>^ = <I)“<I>“. The value of the vacuum term it is not known but, as 
we will see below, this is not really relevant for our purposes. We stress once more that we 
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are not identifying this potential with the SM potential, we just take it as a template for our 
investigations. 

We assume that the metric of spacetime is isotropic and homogeneous, namely 

{t)5ijdx'^dx^ . (2.5) 

This implies that the symmetry of spacetime “overrides” the gauge symmetry, reducing the 
effective degrees of freedom of the gauge field to just one, according to [18, 19] 

^S = 0, A- = f{t)5t. (2.6) 


This choice guarantees that isotropy and homogeneity of space-time are preserved. Similar 
considerations work for static spherically symmetric solutions of Yang-Mills theories [20-22] . 
With the constraint (2.6) we have 




= 6 



(2.7) 


and 






( 2 . 8 ) 


where the dot stands for a derivative with respect to t and <i>^ = 

We now replace these expressions into the Lagrangian (2.1), we work out the equations 
of motion by variation of the fields N, a, /, and and we set N{t) = 1 at the end. The 
Friedmann equations then read 


= 


1 


H = - 


3M2 

1 


3/2 , , 52^2^2 ^ ^ 


2M2 


2/' , 25V" , ^2 , , , 4 pr 


(2.9) 

( 2 . 10 ) 


where H = dhia/dt and matter and radiation densities satisfy respectively the equations 


Pm — Pm 1 Pr — didpj. . 


( 2 . 11 ) 


The system is implemented by the equation for the gauge field degree of freedom 


f + Hf + 


2g^r , 252/^2 


+ 


= 0 , 


( 2 . 12 ) 


and by the Klein-Gordon equations for each component of 


$“ + 3i76“ + ^^-4-^ + A$“(d>2-$2) = 0 , a =1,2,3. (2.13) 

As anticipated in the introduction, and similarly to what happens in the case studied in [10], 
<1>Q cannot be an exact solution of the last equation, unless / = 0 or (7 = 0 or <l>o = 0. In 
principle, however, the condition d >2 = <f>Q can be achieved in the infinite future, if a diverges 
but / does not. As already explained in ref. [10], an eventual late-time dynamics, dominated 
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by a slowly varying effective cosmological constant, cannot be described by the standard 
slow-roll condition (in fact V~^dV /diverges in the large A^-limit) but rather by a “ultra 
slow-roll” regime, where both and can be neglected in eq. (2.13), while the last two terms 
take over and yields a secular variation of the potential. 

This interpretation is further supported by noting that the right hand side of eq. (2.10) 
is negative definite, hence an exact de Sitter-like solution, with H = const, is not allowed. 
However, one can look for asymptotic solutions characterised hy H = const in the infinite 
future, with a —)• oo and the Higgs field settled at the value = <1>q so that its time derivative 
vanishes. 

We now consolidate these qualitative discussions by studying the system of equations 
and its equilibrium points. 


3 Dynamical system analysis 


A very efficient way to study systems of equations as the one in the previous section, especially 
in the context of dark energy, is to map it into a closed system of first-order differential 
equations by defining a new set of dimensionless variables [1]. This technique is also ideal for 
numerical treatment of the equations. We choose the following variables 


V^aM ’ 

V2Ma 


y = 


w 


V2MHa? ’ 
' VSMaH ’ 



(3.1) 


where we introduced the derivative with respect to the e-folding number N = In a, denoted 
by a prime. The conversion between time derivative and A^-derivative is X = HX' for any 
function of time X. 

The deceleration parameter, defined hy q = —1 — H/= —1 — H'/H, indicates 
whether the expansion of the Universe is accelerated or not. When only one source of matter 
is present q is related to the equation of state parameter w by the equation g = (1 -|- 3uj)/2: 
whenever q < 0 (u > —1/3) the expansion is accelerated. In our case, the acceleration is a 
combined effect of the gauge and the Higgs fields, therefore define an effective equation of 
state parameter such that g = (1 -|- 3a;efr)/2. By using the Friedmann equations and the 
variables above we find that 


q= -(1 -h -h — 3v^ + 3z^ — uP') 


where we used the short-hand notation 


(3.2) 


= {zif + {Z2f + {zPl^ , vP = {wif + {W2f + . (3.3) 

With these expressions, eq. (2.9) can be written as the constraint 

11m = ='i-- + y'^ + + r^), (3.4) 

in terms of the matter density Hm- Analogously, we can define the densities of radiation and 
dark energy respectively as 


Qr = P, (3.5) 

llde = Up -\- , (3.6) 
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so that 


1 — ^de • (^■^) 

By differentiating with respect to N each of the variables (3.1), and by using the equations 
(2.11)-(2.13), we find a closed system of eleven first-order differential equations, given by 

I' = 1(1 - lx ), (3.8) 

x' = {q- l)x - l{w^ + 2 y^) , 
y' = y(9 - 1 + 2x0 > 
r' = {q- l)r , 

v' = v{q + 1) + al{wizi + W 2 Z 2 + W 3 Z 3 ) , 

Wi = Wi{q + lx) + V2lyzi , i = 1,2,3, 
z'i = {q - 2)zi - lwi{V2y + av) , i = 1,2,3, 

where a = y/^X/g is a dimensionless constant whose value will be discussed below. 

Before solving numerically the system, it is crucial to find the equilibrium points and 
to study their stability. The system has an infinite number of fixed points that, however, 
correspond to a discrete and finite set of values for the deceleration parameter, namely q = 
2,1,1/2,0, —1. We now examine each family of equilibria. 

q = 2: in this case the set of fixed points is defined by 

z^ = 1, x = 0, y = 0, v = 0, r = 0, 1 = 0, Wi = 0 , i = 1,2,3. 

Physically, this solution corresponds to / = 0 and V = 0, which means that the Higgs field 
takes its vacuum value = <f>Q on the fixed point. By computing the Jacobian we find that 
the eigenvalues are all real but some of them are vanishing thus the equilibrium points are 
non-hyperbolic and numerical methods are necessary to fully assess the stability. In fact, the 
numerical solution (see below) shows that this is an unstable equilibrium that corresponds to 
a Universe dominated by stiff matter (weff = 1) in the remote past. 

q=l: in this case we have two classes of solutions, namely 

(а) = 1 — = {1 — r^)~^ , y = 0, x = 0, Zi = 0, Wi = 0, i = 1,2,3, 

(б) = 1 — x^ — , 'i; = 0, 1 = 0, Zi = 0 , Wi = 0 , z = l,2,3. 

Both classes are non-hyperbolic since some eigenvalues are vanishing. For the class (a) we 
see that the radiation domination cannot be realised entirely with relativistic matter, which 
would require r = 1 and a diverging 1. For the class (b) instead, r = 1 is allowed provided 
X = y = 0, therefore radiation domination can be achieved with relativistic matter only. 
As for the previous case, x = 0 implies that <I>^ = <hQ. The numerical solution shows that 
this equilibrium is not stable, in fact it is “metastable” because the condition q = 1 (i.e. a 
radiation-dominated Universe with cUgfr = 1/3) can last for few e-foldings, for realistic initial 
conditions. 

q = 1/2: here we find the first hyperbolic equilibrium point, which is located at the origin, 
where all variables vanish. The eigenvalues are all real and have both positive and negative 
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sign, thus the equilibrium is unstable. Physically, this solution corresponds to a Universe 
dominated by matter (wgff = 0). By solving the system numerically with the initial conditions 
set arbitrarily close to the equilibrium point, we find that the matter content decreases in 
few e-foldings while dark energy gradually takes over. The radiation content instead does 
not increase since its linearised equation yields r ~ Similarly to the previous case, 

this fixed point can represent the metastable matter domination of the Universe, which ends 
when dark energy becomes dominant. 

q = 0: we have another non-hyperbolic fixed point located at 

vo^ = 1, x = 0, y = 0, ^ = 0, r = 0, 1 = 0, Zi = 0, i = 1,2,3. 

The numerical analysis reveals that it is an unstable point and, in the global evolution, it 
marks the transition between deceleration (cjgff > ~l/3) and acceleration (cUgfr < ~l/3)- 
Note that q vanishes because of the negative sign of w‘^ in the definition (3.2). In typical 
quintessence models, the field does not appear in the equations of motion, only $ does. In 
our model, instead, the Friedmann equations depend also on $ because of the coupling to the 
gauge degree of freedom / via the constant g. Therefore we are forced to define the variables 
Wi (usually absent in quintessence models), which have a non-trivial dynamics and allow for 
the appearance of the critical point at g = 0. In other words, in our model the takeover of 
dark energy (or more precisely, the transition to the accelerated phase) crucially depends on 
the coupling between the Higgs and the gauge field. 

q = — 1: this is the most relevant equilibrium point as it corresponds to an accelerated Uni¬ 
verse (wgff = — 1). It is hyperbolic and located at ^ 

x = 0, y = 0, v = l, r = 0, 1 = 0, Zi = 0, Wi = 0, z = 1,2,3. (3.9) 

All eigenvalues are real and negative except for one so the point is a saddle point. By 
linearising the system, we discover that it is the variable 1{N) that runs away from the fixed 
point. In fact, by solving the linearised equation for y, we find that the scale factor evolves 
as a// ~ In addition, ~ , which indicates that the Higgs field components tend 

to constant values. This is consistent since u —)• 1 nearby the fixed point, so that 

{^l-^l) = 2^MH^. (3.10) 

Finally, the radiation energy density exponentially decays to zero as it can be easily seen 
from the linearised equation for r. Physically, all this means that the Higgs held tends to 
settle asymptotically on some value (which does not necessarily coicide with the vacuum 
expectation value <ho), dragging the Universe towards an accelerated phase, where the Yang- 
Mills held is constant, radiation and matter fade away, while the potential becomes an effective 
cosmological constant. This hnal state corresponds to the only stable equilibrium point 
therefore we conclude that the EYMH system for a homogeneous and hat Universe has a hnal 
dark energy-dominated phase. 

Before proceeding into the discussion of the numerical solutions, we must assign a rea¬ 
sonable value to a, which is in fact the only free parameter of the theory. We note that the 

Strictly speaking, there are two fixed points with v = ±1 but, physically, v is dehnite positive for an 
expanding Universe. 
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stability of the equilibrium points is independent on this parameter since the eigenvalues do 
not depend on a. This can be qualitatively understood from the fact that in the system (3.8) 
a appears only in non-linear terms. Therefore, the structure of the fixed points in unaffected 
by the value of a. However, for numerical solutions we need to estimate it. 

In the Higgs sector of the SM, the coupling constants are related to the vacuum energy 
and to the mass spectrum of the theory. In particular, g = 2Mi4//u, where ~ 80 GeV 
is the mass of the bosons and v = ~ 246 GeV is the vacuum expectation 

value of the Higgs, determined by the Fermi coupling constant. The self-coupling of the Higgs 
is determined experimentally to be A = 0.13. If we identify the parameters g and A in our 
model with these values, we find a = 0.96 With no other available guiding principle, we 
then fix a = 1 in our numerical calculations and we will discuss the physical observables at 
the end. 

4 Numerical solutions 

In the previous section we have shown the existence of one stable asymptotic equilibrium 
point that corresponds to a final state of dark energy domination. In this section we solve the 
dynamical system numerically in the attempt to build a realistic scenario. We anticipate that 
the numerical solutions shown below are just a preliminary study of the full solution space, 
whose exploration requires sophisticated numerical tools and will be presented in another 
work. Our aim is just to demonstrate the validity of our model in recreating a realistic 
scenario. 

We choose to fix the initial conditions at the present time, by imposing Hr ~ 0, ~ 

0.69, and Hm ~ 0.31. In detail, we set 

a: = 10-^ y = 10-8, u = 0.83, r = lO-^, 

zx = -10-^1 , Z 2 = 2x 10-^1 , ^3 = 3 X 10-^1 , ^ = 4 X 10“^ , (4.1) 

wi = 10-9 , ^2 = 4 X 10-1° ,W3 = 2x 10-9 , 

at V = 0. Our numerical solutions show that, with these initial data, matter-dark energy 
equality occurs at N cz. —0.3 (corresponding to a redshift z ~ 0.3) while matter radiation 
equality is at iV ~ — 8.1 (redshift z = 3360). 

The first relevant plot is in Fig. (1) that shows the behaviours of the deceleration param¬ 
eter and of cUgfr as a function of N over several e-foldings. We see that the remote Universe 
starts with a stiff matter-dominated phase (cjgff = 1) and evolves towards the equilibrium 
point at g = — 1, which corresponds to ujeS = ~1- The transition between the two extrema 
is characterised by two plateaux corresponding to radiation and matter domination respec¬ 
tively. They also coincide with two unstable fixed points (<7 = 1 and q = 1/2) discussed in 
the previous section. This proves that the Yang-Mills Higgs equations coupled to gravity 
and ordinary fluids are able to describe all phases of the evolutions of the Universe (except 
for the initial inflationary expansion). Initial conditions may change the duration of each 
intermediate phase but, in general, do not modify the overall evolution of the solutions. 

The evolution of the relative density of radiation, matter, and dark energy is displayed 
in Fig. (2). If we extend the computation at earlier times, we find that the stiff matter 
domination is caused by the variable which becomes dominant in the remote past, see 
Fig. (3). It is also interesting to see the behaviour of other variables. In particular, we show 

^It is amusing how close this value is to the inflationary scalar spectral index Us- 
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in Fig. (4) the plot of zi as a function wi in the present epoch (the results for the other two 
couples of variables are qualitatively the same). The plot clearly shows the attractive nature 
of the fixed point with q = —1. 



Figure 1. Evolution of the deceleration parameter q{N) and the corresponding effective equation 
of state ujes{N)- The initial conditions are given by (4.1) at = 0, corresponding to the black 
vertical line. We note that the evolution of the system mimics the transition between a remote past 
stiff matter-dominated Universe and a hnal dark energy-dominated Universe, passing through two 
transient phases of radiation and matter domination respectively. 
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Figure 2. Evolution of the density parameters with the same initial conditions as in Fig. (1). The 
vertical black line is the present time. 





















































Figure 3. Evolution of the “stiff matter” density parameter + z| + z|, which is dominant 

in the remote past. We find that the time at which it becomes dominant only depends on the initial 
values given to Zi at TV = 0. By lowering these values we can push back the stiff-matter phase at 
virtually arbitrarily large time |7V|. While mathematically interesting, this early time solution is not 
compatible with observations, as it does not connect smoothly to an inflationary phase. 



Figure 4. Evolution of the variable z\ as a function of wi in the range N = [—0.03, 0]. 


As explained above, in principle we do not have clear indications on the value of a. 
However, we can relate it to physical observables. By using the definitions of the variables 
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1{N), y{N), and v{N), we find that 


(4.2) 


M _ y/2 

li 

_ 6V2v 
A'P ayP' 

The first relation implies that, at the present time, and with the initial conditions (4.1), the 
coupling g must be very small. With H = lA x GeV and M = 2.4 x 10^® Gev, and 

a = 1, we find g ~ 10“^^. For a = 1, this implies an even smaller A, of the order A ~ 

The second relation instead gives an estimate of the displacement from the vacuum value of 
the Higgs field, which corresponds to ($^ — $g)/M^ ~ 10 ® with (4.1). 

If we wish to increase the value of g and A, we need to increase of several orders of 
magnitude the initial value of 1{N) (on the contrary, the value of y{N) cannot be larger 
than unity because of the constraints (3.5) and (3.7)). In turn, this implies that the initial 
value of x{N) must be carefully chosen so that the first equation of (3.8) does not yield 
large derivatives. This means, essentially, that x ^ 1/1 over all the integration range. All 
these refinements require more powerful numerical integration methods and will be discussed 
elsewhere. However, the solutions shown above are promising and we believe that a thorough 
analysis will be able to find stable numerical solutions with larger values of A and g. 

5 Concluding remarks 

In this paper, we have studied the full Einstein Yang-Mills Higgs system of equation embedded 
in a flat Robertson-Walker spacetime. The remarkable feature of this system is the set of fixed 
points that coincide with various kinds of matter domination, beginning with an unstable stiff 
matter type in the remote past and ending with a stable final stage of dark energy domination. 
The transition between these two eras is characterised by metastable phases of radiation and 
matter domination. Numerical solutions seem to confirm this evolution and we conclude that 
this model is valid explanation of the current acceleration of the Universe. 

Of course there are many issues to be considered. First of all, a full numerical study of 
the system is necessary, in particular to assess the stability against perturbations. In addition, 
the model can be extended to different kinds of gauge groups and/or potentials. In any case, 
our task was to show that EYMH equations were able to reproduce dark energy. We found 
that not only this is true but also that the entire evolution from radiation domination to 
today is compatible with this model. 

It would be interesting to investigate extensions of this model towards the inflationary 
era. This might be in principle possible recalling that the non-minimal coupling of the Higgs 
field or of the gauge potential to gravity provides for valid models of inflation [14-16]. Typ¬ 
ically, the non-minimal coupling becomes negligible at low energy so it would not affect the 
evolution of the Universe after inflation, which is well described by our model. The possibility 
that the EYMH Lagrangian non-minimally coupled to gravity connects the current accelerate 
expansion to the initial inflationary phase is therefore an open and interesting question. 
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